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Abstract 

The number density (n e ) electron temperature (T e ), 
and velocity of the exhaust of a single pulse thruster 
were experimentally determined for various applied 
axial magnetic fields (0, 1.0, and 2.0 T). At an axial 
station 0.2 m from the anode, the no-field exhaust 
profile is approximately constant up to 4. 0 cm from the 
centerline, n e = 1.0X10 21 p/m 3 , T e =7.2 eV. At 1.0 
and 2.0 T, centerline n e < 10 2 * p/m 3 and the n e 
raises to a maximum of 1.0xl0 21 p/m 3 and 0.'6xl0 2 * 
p/m 3 respectively at a radius of 2 cm. Corresponding 
velocities (0.2 0.8x10^ m/sec) were measured and 

are discussed. 

Introduction 

This report describes the experimental investigation 
of some of the exhaust characteristics of a single-shot, 
pulsed, MPD-ARC thruster. The high power level 
(MW) of the experiment provided an exhaust plasma 
density which was great enough to allow diagnosis using 
Thomson scattering of pulsed laser light. This technique 
gives number density and electron temperature of the 
exhaust without the complication of other auxiliary 
measurements or assumptions. As well, the excellent 
time and space resolution of this technique allow for 
obtaining exhaust profiles (spatial and temporal). 

The exhaust characteristics of the MPD-ARC thruster 
(number density, electron temperature, and velocity) 
reported herein are an essential part of the basic 
knowledge required: 

(1) For feasibility as a thruster, 

(2) For feasibility as a source for a plasma blow- 

down tunnel, and 

(3) For understanding the physics of the exhaust of 

arcs employing auxiliary magnetic fields. 

Pulsed, megawatt level, MPD-ARC investigations 
have been reported by Jahn et al. (•*■) for the self field 
case, the thruster being powered by a delay-line capaci- 
tor bank. The investigation at Lewis Research Center 
of pulsed, megawatt level, MPD-ARC s has been initially 
reported by Michels( 3 , 3). These reports(2,3) describe 
terminal characteristics for the auxiliary magnetic field 
case, the thruster being powered by a crowbarred 
capacitor bank. The arc current decays with time and 
with this type power source. A portion of the detailed 
investigation described by Jahn et al. (*■) gives the 
number density and temperature for one self field case. 
These are obtained with a local probe in direct contact 
with the plasma. Herein, plasma diagnosis through 


Thomson scattering is used as an independent and direct 
method of obtaining number density and electron tem- 
perature for various values of auxiliary magnetic field. 
The technique of determining number density and 
electron temperature from Thomson scattering of the 
beam of a Q-spoiled ruby laser is now well estab- 
lished^* 5 ) . The combination of low electron tempera- 
ture (a few electron volts) and low number density 
(approximately 10 20 particles per cubic meter) have 
been the criterion that determined the simple, yet 
judiciously designed, two-channel, 90° scattering 
system used in this experiment. 

Apparatus 

Capacitor Bank 

The plasma source was energized by a 10 kilojoule 
bank. The equivalent circuit of the bank, crowbar 
switch, and the source is shown in Fig. 1. After the 
bank switch was closed, arc current was allowed to 
develop to its peak current. This takes 21 g see. Then 
the crowbar switch was closed, forcing the arc current 
to decay monotonically with time. The L/R decay 
time ranged from 250 to 350 gsec, depending on arc 
resistance. This allowed an almost linear decay of arc 
current for about 100 gsec after crowbarring and it is 
during this period that the data was gathered. 

Plasma Accelerator 

A photograph of the plasma accelerator is shown in 
Fig. 2. The arc chamber is hidden from view at the 
center of the toroidal dewar for the superconducting 
magnet. The magnet is used to supply the auxiliary 
magnetic field for the accelerator. Plasma flows to the 
left from the accelerator into an evacuated glassware 
system. 

Prior to operating the source, the magnet dewar was 
filled with liquid helium, and after it became cold the 
magnet wire became superconducting. Then the magnet 
was charged to a given magnetic field setting (0 to 2. 2 
T) and maintained at that condition for the testing period. 
The capacitor bank was charged next. It was not 
switched until propellant had properly filled the arc 
chamber (650 gsec). For operational convenience the 
cathode was not preheated for the experiments to be 
described here. 

Nitrogen propellant is introduced into the arc 
chamber by a high-speed gas valve that is actuated by 
an electromagnetic system. All tests were rim at 
0.007 kg/sec nitrogen. The transient, cold flow, gas 
pressure in the arc chamber was measured by a 
piezoelectric pressure pickup in a previous experiment. 
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Using that pressure and orifice equations for steady 
flow, the mass flow rate for all the tests of this report 
were calculated. From the transient pressure records 
it was found that stable flow occurred after 650 psec and 
the arc was started at that time. The arc was formed 
between the cathode and anode ring. A transient plasma 
flows for a few hundred microseconds into the evacuated 
region to the left of Fig. 2. 

A cross-sectional view of the arc chamber is shown 
in Fig. 3. An iron filings map of the magnetic field is 
also shown. The cathode is a tungsten ribbon measuring 
1 cm wide, 2 cm long, and 1 mm thick. The anode is a 
4. 2 cm inside diameter copper ring. 

A sequence controller controls gas puff injection, 
delay for gas distribution, bank switch closure, crowbar 
switch closure, and then data gathering "start" and 
"stop" times. 

In strumentation 

"Thomson Scattering". - Fig. 4 gives the schematic 
of the laser diagnostic equipment and its relationship to 
the plasma test section. The laser was a Q-switehed 
ruby with 2 joules output. The portion of the light 
measured was that which was scattered in a direction 
normal to the plane of polarization of the laser beam 
(i.e. , 90° Thomson scattering). To minimize the stray 
light entering the detection system a series of light 
stops were placed along the path of the incoming beam 
and two light dumps were used. One light dump (no, 1) 
serves as a dump for the main beam, and the other 
(no. 2) serves as a black background for the scattered 
light detection system. The Doppler-shifted scattered 
light was analyzed through the use of a two-channel 
interference filter system similar to that of Patrick. 

The photomultipliers are EMI 9558B's. Each channel 
has a filter holder containing up to seven filters which 
can be easily changed. In the experiment the scattering • 
parameter « £= sin 0/2)] is much less than 

one («* 0 . 1). A is the laser wavelength and A^ the 
Debye distance. Theoretically then, if a Maxwellian 
electron velocity distribution is assumed, single-shot 
data from two filters with different peak wavelengths is 
sufficient to determine the temperature. Using the 
available transmission curves for the interference 
filters the fraction of scattered light passed by each 
filter was calculated as a function of electron tempera- 
ture. The ratio of these fractions for filters with peak 
wavelengths of 695. 0 nm and 697. 1 nm is given in Fig. 5. 
In taking data the ratios were computed for each shot and 
then a series of five ratios was averaged before deter- 
mining the temperature from Fig. 5. 

Temperatures could also be determined by using just 
one of the channels. In this case ten shot averages were 
taken for a series of six filters ranging in peak wave- 
length from 694. 3 nm to 697 . 8 nm. A best fit was then 
made of the theoretical curves to determine the temper- 
ature. Fig. 6 shows a comparison of theoretical and a 
typical experimental profile. It was found that temper- 
atures obtained through the two methods agreed within 
the limits set by nonreproducibility of the discharge and 
other statistical factors. 


Most of the data herein was taken using the two- 
channel method, since fewer shots were needed. There 
was also the added advantage that signal ratios taken 
from the same laser shot and plasma discharge were 
more accurate than tho se derived from two separate 
shots. 

The system was calibrated for density by measuring 
the Rayleigh scattered light when the test section was 
filled with nitrogen at a few millimeters pressure. The 
shot to shot variation in the calibration signals was less 
than ±10 percent. 

For scattering from the plasma, shot to shot varia- 
tions ranged from ±10 to ±30 percent. Examination of 
other diagnostics such as light emitted over the plasma 
and Faraday cup signals indicates that this shot to shot 
variation in Thomson scattering can be explained by the 
nonreproducibility of the discharge. 

The lower limit in density which could be accurately 
determined in this experiment was about 0. 5 X 10^0 m“3. 
This limit occurs because of the noise on the photo- 
multiplier signal resulting from the light emitted by the 
plasma. This limit can obviously be reduced by the use 
of a higher powered laser. Fig. 7 shows typical two- 
channel scattered light signals. 

"Faraday Cup Probe" . - A Faraday cup probe was 
used to collect ions from the streaming plasma simul- 
taneous with scattered light data gathering. It was 
biased with a 45 V battery to repel electrons from the 
cup. The probe face was kept 10 cm downstream from 
the scattering measuring station so as not to introduce 
stray light signals in the scattering chamber. The 
probe could be positioned to be on the same duct radius 
as the scattering measurement. Pinholes could be 
changed on the probe face (0.0002 m dia. or 0.00074 m 
dia. ). The probe system was kept electrically floating. 
Instrument isolation amplifiers with excellent common 
mode rejection were employed. Si order to obtain iso- 
lation and good common mode rejection instrument 
response time was compromised at 15 p sec rise time. 
Thus the leading edge and fine detail of the cup signal 
were not available in this experiment. Since collected 
ion current is proportional to area, a check was made to 
insure that this was the case by changing hole size. This 
proved to be linear. 

Results and Discussion 

The value of a direct method of determining number 
density through Thomson scattering was demonstrated 
clearly in the early stages of this experiment. It was 
found that minor changes in the insulator design near the 
cathode in the chamber (all other dimensions and mate- 
rials remaining the same) seriously affected the exhaust 
number density. Factors of two to three increase in 
number density were obtained by fully insulating the 
cathode supports. This forced the discharge to be 
initiated and seated at the cathode ribbon. In earlier 
designs, a very minor arcing damage was noted on the 
cathode supports. This probably occurred only briefly 
in random arc initiations. In this report, only the data 
for the most improved version will be described. 
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The precise temporal and spatial resolution of the 
Thomson scattering technique identified another feature 
of the exhaust. Slot to shot variations in number density 
and temperature were serious. The variations in the 
manner of initiation and development of the arc thus were 
observable as exhaust variations. The data presented is 
the mean value of at least five shots for each condition. 

The exhaust flowed from the source into a 15 cm 
diameter evacuated glassware duct. A detailed examin- 
ation of the exhaust was undertaken at a station 20 cm 
downstream from the anode face. The number density 
and electron temperature profiles were measured at the 
station for various auxiliary magnetic field cases (0, 

1.0, and 2.0 T). The results are shown in Fig. 8 for 
the case where peak arc-current for each shot is 13.4 
kA. All shots had the same mass flow rate (0. 007 kg/sec 
N 2 ), and the same gas settling time of 650 microseconds. 
The data was taken about 20 microseconds after the 
plasma light front had reached the measuring station. 

This time was chosen so that initial transient effects 
would be small and the measurement would be made in 
the streaming plasma portion of the flow. 

The time history of the plasma light passing the 
measuring station has a general shape shown below in 
Fig. 9. Shot to shot this shape varies, sometimes the 
pedestal portion is distinct, at other times it is part of 
the smoothly decaying tailoff portion. Generally, the 
spike portion has less number density by at least an 
order of magntiude. The pedestal has N e values in the 
range of 10 20 - 10 21 p/m 3 and the tail portion has 
lower and time-decaying values. Moving the measure- 
ment time later (to 60 - 7 Op sec on the sketch) quite often 
showed factors of two reduction in number density. This 
is to be expected since during that corresponding time 
the arc power is reducing also by like factors due to the 
decay in arc current. 

With these general characteristics understood, the 
data of Fig. 8 can now be discussed. For the case of no 
auxiliary field, the number density is 1.0x10 p/m . 

The density profile was nearly constant from duct 
centerline to 4 cm radius. The electron temperature 
falls linearly from a centerline value of 7. 3 eV to 3. 5 eV 
(4 cm). The arc power was 2.3 MW at 36 psec after the 
arc was initiated. That is 15 psecs after the crowbar 
time. The number density profile was much different 
for the cases where an auxiliary magnetic field is 
employed. Over an order of magnitude reduction in 
density occurs on centerline for both the 1.0 and 2.0 T 
cases. The peak density occurs at a radius of 2. 0 cm. 

At the 2 cm radius for the 1.0 T case, the number 
density was still about 1.0X10 2 -*- p/m 3 but for the 2.0 T 
case it reduced to 0. 6x10 21 p/m 3 . For the 2.0 T case, 
the arc power was 9. 5 MW, at 36 psec after the arc was 
initiated, over four times greater than the "no field" 
case. This power increase was not manifest by number 
density increase, nor electron temperature increases 
from the "no field" case, in fact, the density reduced! 
With these facts, the next logical exhaust characteristic, 
velocity, was measured. 

The time of arrival of the light front at the measuring 
station was used to determine the velocity, Vlite> 


shown in Fig. 8. This is an indicator of light front 
behavior only and not necessarily ion velocity after the 
front has passed. For the "no field" and 1. 0 T cases, 
Vute were alike, within measurement error at about 
0.4x10^ m/sec. The velocity for the 2.0 T case was 10 
percent greater. This would not be a great enough in- 
crease for the power increase to be explained as in- 
creased streaming energy. 

An attempt was made to obtain ion velocity through the 
use of a biased Faraday cup. At the 13. 4 kA peak cur- 
rent case, only one condition was examined success- 
fully. This was because at higher power levels the ratio 
of the common mode signal (to the measurement signal) 
could not be reduced enough. From the one successful 
cup measurement at 13.4 kA and the corresponding 
number density, the ion velocity, Vj, was found to be 
0. 58x104: m/sec. This prompted investigation of a lower 
power level regime where the cup could be used succes- 
sfully. 

At the lower power regime (7.4 kA peak arc current) 
the measurements were made at r = 2cm (where N e 
is the greatest value). The data are shown for the 
7.4 kA peak current case in Fig. 10 and the 13. 4 kA data 
is included for comparison. The number densities were 
lower for the 7. 4 kA case, as are the front velocities. 

Ion velocities, from the Faraday cup measurements, 
showed a six-fold increase as the field is increased from 
0 to 2 T. Thus, at this power level, the ion velocity 
increase comes closer to being proportionate to the 
observed power increase. Again, as for the 13.4 kA 
cases, it was found that Vjjjig increased slightly with 
applied magnetic field. Thus Vute is only an indicator 
of front arrival effects and not the streaming ion flow 
component of the exhaust. 

The product N e Vj should indicate the exhaust 
momentum for that station and time. This product, when 
plotted versus applied magnetic field (see Fig. 11) peaks 
at about 1. 0 T for the data of Fig. 10. 

The product N e Vj 2 should be proportional to dynamic 
stream pressure. For the data of Fig. 10, N e Vj 2 is 
almo st linearly increasing with applied magnetic field. 

It is planned, in the near future, to examine, with an 
independent measurement of dynamic pressure, if the 
pressures obtained agree with the Faraday cup derived 
pressure calculations. 

Concluding Remarks 

The number density and electron temperature of the 
exhaust of a pulsed megawatt MPD-ARC plasma thruster 
were determined for two different peak arc-current 
cases (13.4 kA and 7.4 kA) for various values of auxili- 
ary magnetic field (0. to 2.0 T). The radial profile of 
the exhaust was obtained for a position 20 cm down- 
stream from the anode face. For the 13.4 kA case, with 
no auxiliary field, the centerline number density was 
1X10 21 p/m 3 and was relatively flat out to 4. 0 cm 
radius. For the applied field cases the centerline 
number density is reduced to below 10 2 ® p/m 3 and 
peaks at 0.02m radius. At the 2 cm radius location, 
for the 1.0 T case the number density is still about 
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l.OxlO 21 p/m 2 but as field is increased above this point 
number density reduces. At 2.0 T the number density is 
down to 0.6x10 21 p/m3. An electron temperature range 
from 7. 3 to 3.0 eV was noted. . Using the Thomson 
scattering technique is a very powerful aid in examining 
high density exhaust. It gives local N e and T e directly 
without auxiliary methods or assumptions and with 
excellent temporal and spatial resolution without 
perturbing the plasma. 

The same degree of confidence cannot be applied to the 
velocity measurements discussed herein. These suffer 
the usual "probe-in-plasma" uncertainties. The Faraday 
cup derived velocities determined for the 7. 4 kA range 
from 0.1 to 0. 68xlo4 m/sec. The probe perturbs the 
plasma an unknown degree, has the attendant particle 
gathering problems, and cannot be used in very high 
power, extended pulse duration, exhausts. 

The temperature measurements for the "no-field" and 
field cases show that increasing arc power through the 
aid of a passive externally applied, magnetic field is not 
manifest by increased electron heating. Although not 
conclusive, but rather useful to suggesting the direction 
of further investigation, the lower powered part of this 
experiment has suggested that the increased power from 
the "no field" case tofield cases does manifest itself 
primarily in increased plasma streaming energy. More 
work needs to be done to verify this through another 
independent velocity measurement, over more stations 
and at higher power. 
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Figure 2. - Plasma accelerating system. 
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Figure 3. - Arc chamber. 
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Figure 4. - Schematic of scattering diagnostic equipment. 
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Figure 5. - Temperature vs signal ratio. 



FILTER CENTER WAVELENGTH, nm 


Figure 6. - Theoretical and experimental spectral 
profiles. T e ~ 6 eV, n e ~ lxlO 21 m" 3 . 
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Figure 7. - Typical two-channel photomultiplier signals. 
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Figure 8. - Exhaust characteristics. 
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Figure 9. - Typical plasma light time history. 



Figure 10. - Exhaust characteristics vs magnetic field 
at r =0.02 m, z =0.2 m. 
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Figure 11. - Exhaust N e Vj and 
N g V? vs magnetic field strength at 

r = 2 cm, z = 20 cm, 7. 4 kA. 
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